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We review the method of differential renormalization, paying special attention to a new 
constrained version for symmetric theories. 

1 Introduction 

Popular regularization and renormalization methods work in momentum space. Typ- 
ically, the divergences which appear in loop integrals at large internal momenta (ul- 
traviolet divergences) are first regulated — i.e., the integrals are modified so that they 
are finite, but diverge in the limit of no regulator — and then substracted by adding the 
necessary counterterms. For example, in dimensional regularization Q the regulated 
integrals are defined in n dimensions by analytical continuation, with n an arbitrary 
complex number. Eventually, the poles appearing at n = 4 are cancelled by appropriate 
counterterms and a finite (renormalized) result is obtained for n — > 4. Renormaliza- 
tion without intermediate regularization is also possible, as in the BPHZ method Q, 
where the first terms of the Taylor expansion in external momenta of the integrand 
are substracted off before integrating. Although momentum space is more natural for 
calculations of scattering amplitudes with fixed external momenta, nothing prevents us 
from working in coordinate space and, if required, perform a Fourier transform at the 
end. In coordinate space, ultraviolet divergences correspond to a singular behaviour at 
short distances. 

In Ref. H a method of renormalization in coordinate space was proposed: differ- 
ential renormalization (DR). It is based on the observation that primitively divergent 
Feynman graphs are well defined in coordinate space for non-coincident points, but 
too singular at coincident points to allow for Fourier transform. In other words, the 
corresponding expressions are not well-behaved distributions. The idea of DR is to 
substitute the singular expressions by derivatives of well-behaved distributions, in such 
a way that the former ('bare') and latter ('renormalized') expressions are equal almost 

1 Presented at the XXI International School of Theoretical Physics "Recent progress in theory and 
phenomenology of fundamental interactions" , Ustroh, Poland, September 19-24, 1997. 



Figure 1: One-loop diagram contributing to the four-point vertex in <£ 4 . 



everywhere. These derivatives are prescribed to act formally by parts in integrals such 
as Fourier transforms. In this way, finite Green functions are obtained, without the 
need of intermediate regularization. DR acts directly on bare Feynman graphs and 
does not introduce explicit counterterms. The procedure is best illustrated in terms of 
one example: the one-loop four-point bubble graph of massless <3? 4 (Fig. [I]). We work in 
euclidean space, which leads to simpler functions. The massless propagator in position 
space is A(x — y) = ( x l y yi and the vertex, — \5(x\ — X4)5{x2 — x^)5(xz — X4). The 
bare expression for the amputated graph is 

A 2 1 1 

r(x 1 ,x 2 ,X3,x±) = YlQ^ 5 ( Xl ~ X2 ^( X3 ~ x ^ fa _ X3 )4 + 2 P erms - U) 

This involves the singular function which has a logarithmically divergent Fourier 
transform. To renormalize it with DR, one must solve a differential equation and find 
f(x 2 ) such that 

^ = °/(* 2 ) (2) 

for x 7^ 0. Actually, one derivative would be enough in this case, but one uses the 
D'Alambertian □ = d^d^ to preserve manifest euclidean invariance. The solution of 
Eq. ©is 

2 l \ogx 2 M 2 

f( X ) = "4 ~2 > ( 3 ) 

where M is an arbitrary constant with dimensions of mass, required for dimensional 
reasons, and we have omitted a possible but irrelevant additive constant. Although we 
shall not discuss it here, it is worth mentioning that the constant M plays a central role 
in DR: the renormalized amplitudes satisfy renormalization group equations, with M 
the renormalization scale. The renormalized expression of the singular function reads 
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With the formal integration by parts rule, this is a tempered distribution which admits 
a finite Fourier transform: 
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where M = 2M/je, and je = 1.781... is Euler's constant. Substituting Eq. 
Eq. (||), the renormalized vertex graph is obtained: 
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The renormalized expression in momentum space follows from Eq. (B): 
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DR has been successfully applied in different contexts: the Wess-Zumino model Q, 
lower-dimensional || and non-abelian gauge theories g, two-loop QED @, a chiral 
model ||, a non-relativistic anyon model |J, curved space-time and finite tempera- 



ture [ 10 1 , the calculation of (g — 2)i in supergravity |11], Chern-Simons theories [12] 



and non-perturbative calculations in supersymmetric gauge theories [13]. Other for- 
mal aspects of the method have been developed in Refs. [^, [R], 17] and different 
versions of DR can be found in [18, [ig|] . 

When symmetries are an issue, it is important that the renormalization program 
preserves the corresponding Ward identities. In general, even when the regularization 
procedure breaks some relevant symmetry, one can still recover it by the addition of 
finite counterterms^. In practical calculations and formal proofs to all orders, it is 
nevertheless more convenient to have a method that directly preserves the Ward iden- 
tities. The great success of dimensional regularization is mainly due to the fact that 
it automatically respects gauge invariance. It is known, however, that it has problems 
in dimension-dependent theories like chiral and supersymmetric theories. Its variant 
dimensional reduction |2(J is usually employed in these cases, although inconsistencies 
may arise at high orders [21]. DR does not change the space-time dimension and it was 
expected to become a renormalization procedure respecting gauge and chiral symme- 
try. In its original version, however, this is not automatic. The ambiguities inherent 
to the manipulation of singular functions are taken care of by introducing arbitrary 
renormalization scales for different diagrams. Different choices of the renormalization 
scales give rise to different renormalization schemes and only a subset of these schemes 
corresponds to a symmetric renormalized theory. Hence, the scales must be fixed to 
enforce the relevant Ward identities. A change of renormalization scales is equivalent 
to the addition of finite counterterms, so the situation does not differ much from the 
one with symmetry-breaking regulators. 

In Refs. [O, 22 1 a procedure of DR was proposed which fixes all the manipulations 



and only introduces the necessary renormalization group scale. This constrained DR 



2 The exception is called an anomaly: the quantum renormalized theory does not have a symmetry 
of the classical theory. 
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has been explicitly shown to respect the one-loop Ward identities of abelian gauge 
symmetry [22|fj and to preserve supersymmetry in a supergravity calculation [11]. 

In the following, we first describe the method of constrained DR and then use it 
to calculate in detail the electron self-energy and the vertex correction in QED, as an 
illustration. We also derive the corresponding momentum space expressions and check 
that the corresponding Ward identity is automatically fulfilled. 



2 Constrained differential renormalization 

In this section we briefly describe the constrained procedure of DR at one loop intro- 
duced in Ref. . The idea is to find a consistent way of performing the manipulations 
of singular expressions carried out in the process of renormalization. It turns out that a 
small set of formal rules is sufficient to completely fix the renormalization scales (except 
one, associated with the renormalization group invariance). Furthermore, the resulting 
renormalized amplitudes were explicitly shown in Refs. |22| , |TTf to satisfy the one-loop 
Ward identities of abelian gauge symmetry and to render a vanishing value for the 
magnetic moment of a charged lepton in supergravity (which is the required value if 
supersymmetry is respected p4|). 

The set of rules contains the two basic DR rules: the use of DR identities like Eq. (Q) 
(always with the same renormalization scale!) (rule 1) and the formal integration by 
parts prescription (rule 2). In addition, we need another two rules. One is (rule 3): 

[F(x,x 1 ,...,x n )5(x - y)] R = [F(x,x 1 ,...,x n )] R 5(x - y) , (8) 

where F is an arbitrary function. The other one requires the general validity of the 
propagator equation (rule 4)'- 

F(x,xi,...,x n )nA(x) = F(x,x 1 ,...,x n )(-5(x)) , (9) 

where A(x) = -^x-^i is the massless propagator (the massive case is analogous). This 
is a valid mathematical identity between tempered distributions if F is well-behaved 
enough. This rule formally extends its range of applicability to an arbitrary function. 
The main point of the constrained method is to require consistency of renormalization 
with these rules. Such requirement fixes all the ambiguities in DR at least to one loop. 
Let us explain how the rules are used in practice with some simple examples. The 
results will be used in the next section. We first introduce some convenient notation: 
we define the bubble and triangular basic functions as 

B[0] = A(x)OA(x) , (10) 
T[0] = A(x)A(y)O x A(x-y) , (11) 

where O is a differential operator. The significance of this kind of functions is that 
any one-loop bubble or triangular Feynman diagram can be expressed in terms of them 

3 In Ref. Smirnov presented an abelian gauge invariant method within his version of DR. 
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(and their derivatives) using only algebraic manipulations and the Leibnitz rule for 
derivatives. Hence, the problem of renormalization reduces (at this order) to finding 
the renormalized expressions of the singular basic functions. Note that B[0] is singular^] 
at x = for any 0, whereas T[0] is only singular (at x = y = 0) when O contains 
two or more derivatives. These basic functions are easily renormalized using the set of 
rules described above. For example, 
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For basic functions with non-trivial tensor structure the procedure is more involved 
and can be found in Ref. In particular, we shall need in the following the identity: 



+ 
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where the first term is finite thanks to its tracelessness. The local term proportional to 
6(x)5(y) was not considered in the earlier literature (before Ref. p2| , where Eq. ( |i~4| ) 
is worked out in detail) and comes from imposing consistency with the propagator 
equation (rule 4). Notice that 



(15) 



This might seem strange, but in fact also occurs in other schemes like dimensional 
regularization or Pauli-Villars. Differentially renormalized expressions of basic func- 
tions appearing in one-, two- and three-point one-loop Green functions can be found 



in the tables of Ref. [22|. The treatment of four-point one- loop Green functions will be 
presented in Ref. [25]. 



4 In this paper the term 'singular' should always be undertood as 'too singular to allow for a Fourier 
transform'. 
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Figure 2: Feynman rules of massless QED in euclidean coordinate space. 

x 3 




electron self-energy vertex correction 



Figure 3: One loop diagrams contributing to the electron self-energy and the electron- 
electron-photon vertex in QED 



3 Simple applications 

We now show how the method works in practice with two detailed examples: the 
renormalization of the one-loop 1PI electron self-energy and electron-electron-photon 
vertex in massless QED. The Feynman rules of massless QED in euclidean coordinate 
space are gathered in Fig. with {7^,7^} = 25^. We work in the Feynman gauge. 
The case of an arbitrary Lorentz gauge was discussed in Ref. p^]. Let us first calculate 
the electron self-energy, given by the first Feynman graph in Fig ||. The bare expression 
is 

E(x) = e 2 7Q A(x) ^A(x) 7a , (16) 

where x = x\ — X2- Due to translation invariance, £ only depends on x. Notice that 
Eq. ( |l6| ) involves no integration, in contradistinction with the corresponding expression 
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in momentum space^J. After some straightforward (four-dimensional) diracology and 
the use of Leibnitz rule to extract the derivative, one obtains 

S(x) = -e 2 ^(A(x)) 2 

= -e 2 0B[1] . (17) 



The renormalized value is, from Eq. (12), 

p . 1 9 „ logx 2 M 2 , . 

Let us now deal with the vertex correction (see Fig. ^). Reading directly from the 
Feynman rules, 

V„(x, y) = (-ie) 3 7a fA^f- ^)A(y) la A(x - y) , (19) 

with x = xi — X3 and y = X2 — 23. Simplifying the Dirac algebra and using systemat- 
ically the Leibnitz rule to rearrange derivatives, V^(x,y) can be expressed in terms of 
triangular basic functions: 

V,(x, y) = ie 3 {-2 7fe7M7a (^^T[l] + d x a T[d b ] - dfT[d a }) 

-2 7M T[D] + 4 7a T[3 a d M ]} . (20) 

The renormalized expression is obtained directly from Eqs. (|l3|) and (|i~4]): 



V^foy) = ie 3 {-2 7Wa (^^T[l]+^T[a 6 ] 
- ^T[0 a ]) + 4 7a T[d a 3 M - - A 5 a ^u) 

1 1 logx 2 M 2 . 11 ... 

" 4 (4^2)2 >° ^2 " ^ " s^^^)^)} • ( 21 ) 

Once the graphs have been renormalized in coordinate space, one can perform a Fourier 
transform (without any regulator) to obtain the corresponding finite expressions in 
momentum space. We need the Fourier transforms of the basic functions in Eqs. (18) 
and (^T|). The latter are more involved, so let us see in some detail how to calculate 
them. The Fourier transform of a distribution "of two variables", /(x,y), is 

/W) = / d 4 xe l ^e i yP'f(x,y) . (22) 

With the integration by parts prescription, total derivatives in f(x, y) yield 

f%-> -ip, ; 8* - -ip', ■ (23) 

5 In general, coordinate space calculations involve one integral less than the momentum space ones. If 
one is interested in quantities defined at determined momenta (like scattering amplitudes), the missing 
integral must be done at the end as a Fourier transform (without regulator!). One example when this 
Fourier transform is not needed is the calculation of beta functions. 
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For the finite triangular functions we have: 

T[0] = J d 4 xd\je ix - p e iy p 'A(x)A(y)O x A(x-y) 

d i xd%e lxp e iy - p ' 

d% d% d 4 k e _ ix . kle _ iy . k2e _ i(x _ y) . k _ 6{k) 



(2vr) 4 (2vr) 4 (2vr) 4 k\k 2 k 2 

d% d% d A k 6{k) 



(2^) 4 (2vr) 4 (2vr) 4 k\ k 2 k 2 
x (2ir) 4 5(p - Jfei - k)(2n) 4 5(p' -k 2 + k) 

d A k 6(k) 

(2T) 4 (p- k) 2 (p' + k) 2 k 2 ' 



(24) 



where O(fe) is obtained from O 31 by the replacement d x — ► —ik. The integrals in 
Eq. ( pip appear (with a regularization which is not present here) in standard one-loop 
calculations in momentum space and can be evaluated with standard techniques. We 
shall do that later on, in the limit p' — > 0. On the other hand, the Fourier transform of 
the renormalized basic function in Eq. (|2l] ) reduces to Eq. (||): 

!■*[□] = / (&dfr e^jvP' -L □ bg X l M2 Six - y) 

J 647T 4 X A 

= [d*xe-< p ^'^ x2M2 
J 64tt 4 



x 2 



1 (25) 



16vr 2 faV M 2 

With all these formulae, we get the renormalized vertex correction in momentum space: 



V*(p,p') = ie 3 {- 7bl ^ a (- Pa p' b f[l] - i Pa t[d b ] + ip' b f[d a ]) 
+ 4 7a f[a o ^-^ aAt D] 



l 67r 2 7,(log(^^) + 2)}- (26) 

The Fourier transform (in one variable) of the self-energy in Eq. (|l8|) is directly given 
by Eq. (§: 

£*(p) = -i^ ^log(^) . (27) 

4 The vertex Ward identity 

Finally, let us verify the Ward identity between the renormalized vertex correction and 
electron self-energy. For 1PI Green functions it reads 

(d x + d»)V*{x, y) = zeS R (x - y)(6(x) - %)) , (28) 
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where d X3 f(xi — x^,X2 — 23) = — (d x + d y )f(x,y) has been used to express it in the 
translated variables x and y. At points away from the origin this identity must hold, 
since the bare Green functions have not been modified there. A possible disagreement 
can only arise from the contact terms at x = y = 0. In fact, both sides of Eq. ( p8|) 
are distributions, and to compare them one must either use formal properties of delta 
functions, etc, or integrate with an arbitrary test function 4>(x, y). In particular, one can 
perform a Fourier transform (<p(x,y) = e lx ' p e iy ' p ), which retains all the information. 
In other words, we can check the Ward identity in momentum space, 

- i(p, +P',)% R (P,P) = ie[£ R (-p') - £ R (p)] , (29) 

using the momentum space renormalized Green functions in Eqs. ( |2"6| ) and (|27|). For 
simplicity, we consider the limit p' — * (i.e., the Fourier transform in y reduces to an 
integral without any weight). In this limit the relevant integrals reduce to 

*»J -° IS?F ' (30> 
t M - '3> - JL. ™>i -J^l , (31) 

whereas T[l] is logarithmically infrared divergent and p^T[l] — > 0. With these values, 
we obtain for both sides of the Ward identity in Eq. fl29|), in the limit p 1 — > 0, the same 
result: 

Since both sides are equal, the Ward identity is indeed satisfied: constrained DR has 
preserved it automatically, i.e., without any a posteriori adjustment. 

5 Conclusions 

DR is a renormalization method which works in coordinate space and does not introduce 
any intermediate regulator ||]. We have shown how it can be easily applied to the 
calculation and renormalization of one-loop Feynman diagrams. In many cases, worked 



out calculations at higher orders are also simpler in DR than in other methods y, 15]. 
Some nice features of the method are the following: 

• It is minimal, in the sense that the Green functions are never modified except at 
the singular points. 

• It does not change the space-time dimension, easing the diracology and tensor 
manipulation, and the treatment of chiral theories [||, ||, 11]. 



One integration less has to be performed, unless one is interested in some quantity 
defined for fixed external momenta. In such case, one has to Fourier transform 
the renormalized expressions, but without any regularization. 
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Some overlapping divergencies disentangle in coordinate space. 

It is better suited for theories which are naturally defined in coordinate space, 
like theories with conformal invariance Q, in curved space or at finite tempera- 



ture uj 



In the constrained procedure of DR (2^] all the local terms are fixed. This determines a 
renormalization scheme which turns out to be symmetric in all known examples (2^, [llj . 
Here, we have used constrained DR to calculate two one-loop Green functions in QED 
and we have verified that the Ward identity relating them is automatically fulfilled after 
renormalization. We have dealt with massless theories for simplicity. The treatment of 
massive theories in (constrained) DR is worked out in Ref. 17] (Ref. 
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